www.gusucode.com > 基于粒子滤波的故障检测,使用似然函数作为检测函数 > 基于粒子滤波的故障检测,使用似然函数作为检测函数/code/FDI based on SIR likelihood/drawpar.m
function xnext = drawpar(xprev, stdw, k) M = length(xprev); b0=25; for i=1:M %if i<301 % b=b0; %else % b=b0*0.5; %end xstate = xprev(i)/2 + b0*xprev(i)/(1+xprev(i)^2) +8*cos(1.2*(k-1)); xnext(i,1) = stdw*randn+xstate; end